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Recently, Gibelli published a paper [1| about the derivation of the first- and second- 
order velocity slip coefficients. In this paper, the author chooses of comparing his results 
with those obtained by ,„ea„s of two vanationa. methods Q. g underlying how h. 
method predicts results close enough to those reported in [2|], while there is a "remarkable 
disagreement" with those presented in |J|, although the latter "compare favorably with 
the experimental data" . The aim of this comment is to try to trace back the origin o " 
this discrepancy considering that the two variational principles presented in [2] and in 3| 
are indeed deeply different. The variational method used in [2'| applies to the integral 
form of the linearized Boltzmann equation, which is available explicitly only for simplified 
kinetic models. In the more general case of the Boltzmann equation based on the true 
linearized collision operator and general boundary conditions, it is not possible to obtain 
closed form expressions for all the different operators appearing in the integral represen- 
;ation, which should consequently be approximated by series expansions. Therefore, in 
21], the authors obtain approximate solutions of an approximate integral equation where 
truncated series need to be managed in order to obtain solutions of the Boltzmann equa- 
tion in closed form. The last point (that is, the use of truncated series) is the common 
link with the method of solution presented in l|, although within a completely different 
framework. On the contrary, the variational technique presented in |3| applies directly to 
the integrodifferential form of the linearized Boltzmann equation and can be used for any 
linearized collision term and extremely general boundary conditions. This means that. 
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2|], the right Boltzmann equation for hard-sphere molecules has been solved without 
approximations in its form. Only two eight-fold integrals have been evaluated numerically 
by using a Monte Carlo integration, but since these two integrals are of the same simple 
form of those analytically solved in l^^, they could be computed with high accuracy. All 
the rest has been carried out analytically. The method of solution presented by Gibelli 
requires the numerical evaluation of several eight-fold integrals but in this case the sta- 
tistical error in the Monte Carlo integration is not negligible, as himself states (pag. 6 
in [l|). Concerning the role of the test function, there exist several basic theorems which 
allow to perform a good choice. In our context, then, this choice is even simpler since 



the functions which need to be approximated by trial functions are solutions of physically 
realistic problems and expecially in their asymptotic forms can be immediately obtained 
via the use of the Chapman- Enskog procedure, as done in [^ , or via the solution of the 
Boltzmann equation in integral form based on a simplified kinetic model, as done in |3|. 
In l|, pag. 11, the author writes:"... the current approach ... provides the velocity slip 
coefficients with good accuracy in the entire range of the accommodation coefficient" . It 
is difficult to imagine how he could deduce this good accuracy ... The only explanation 
could be that he infers the accuracy of his asymptotic near-continuum solution starting 
from the accuracy of his solution for the Poiseuille mass flux valid in the whole range of 
the Knudsen numbers. Unfortunately, this inference is completely wrong from the math- 
ematical point of view. This is proved by Gibelli himself. Indeed, Fig. 1 and Table II in 
l| show that the results concerning the complete velocity and flux profiles obtained by 
Gibelli are in fair agreement with those presented in |5| for different Knudsen numbers 
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and accommodations coefficients. But then Fig. 3 in [1[ shows clearly that the asymptotic 
solution for the Poiseuille mass flux truncated at second order is completely wrong in the 
slip flow region (where the second-order slip coefficient should be deduced) . A further ex- 
ample in this respect is offered by the BGK model. The solution of the Poiseuille flow rate 
problem obtained by using the linearized Boltzmann equation based on the BGK model 
slightly deviates (only a few percent) from that obtained by using the linearized Boltz- 



mann equation for hard-sphere molecules 



6| , while the asymptotic second-order solutions 



differ considerably from each other, as pointed out in [3| . This is due to the fact that the 
BGK model is a 'flrst order' model. Therefore, it can not be considered so accurate when 
higher order terms are concerned, since a single relaxation time model can not give a good 
approximation for both flrst order and second order effects 7|]. The method of solution 
presented in pj] suffers from the same problem, although within a completely different 
framework. In fact, even if Gibelli starts from the exact hard-sphere Boltzmann equation, 
the method of solution that he proposes consists in expanding the distribution function in 
terms of half-range Hermite polynomials and then a system of moment equations in the 
expansion coefficients is solved so as to "satisfy approximately the hard-sphere Boltzmann 
equation", as himself states. Therefore, in order to flnd closed form solutions, he has to 
truncate the polynomial expansion of the distribution function at some flnite order. The 



final result of the procedure presented by Gibelli in l| is summarized in Fig. 3 of his 
paper, where it is clearly shown that the second-order description of the flow rate Q is 
completely incorrect in the slip region. To overcome this difficulty Gibelli suggests to 
introduce a third slip coefficient, A3, and presents the following expansion for the volume 
flow rate: 

(where the bar is used to distinguish the velocity slip coefficients which account for the 
structure of the Knudsen layer). Then he writes: "the volume flow rate given by Eq. (27) 
with Aq = 1, in fact, can also be obtained by integrating the velocity proflle satisfying 
the linearized Navier-Stokes equation, subject to the boundary conditions" at the walls: 

e. = TA.k/-^ - K^\M - AsK^)^ (28) 

ay ay"^ 

where y is the coordinate normal to the walls, C,x is the tangential velocity component and 

Kn is the Knudsen number. Unfortunately, this inference is incorrect in several respects. 

• From the mathematical point of view, since the slip coefficients which multiply each 
normal derivative of the velocity proflle in the boundary condition (28) should have 
an order of accuracy in the Knudsen number as the order of the derivative itself, that 
is A2 is a second-order slip coefficient and should multiply the second derivative of 
^x, while A^ is a third order slip coefficient and should multiply the third derivative 
of $,x- This is easily deduced by dimensional reasons and is clearly proved in |7|, 

8|] . But in Eq. (28) in front of the second derivative of ^x appears a strange mix of 
coefficients of different order in the Knudsen number. 

• From a physical point of view, since the third derivative of the Poiseuille velocity 
proflle is zero, therefore third order effects in the flow rate Q (included in ^43) are 



not expected to arise from the third order slip coefficient (which is t 
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front of j-(-j%)) but from the contribution of the Knudsen layer [7| 
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From a conceptual point of view, since Gibelli writes the boundary conditions to be 
associated to the Navier-Stokes equation (Eq. (28)) in terms of the 'effective' slip 



coefficients (that is, slip coefficients which account for the structure of the Knudsen 
layer). Instead, a successful higher-order slip model is one that does not agree with 
the Boltzmann equation solutions in the Knudsen layers, in order to be able to 
predict within the same framework all the physical quantities of interest: not only 
the flow rate but also, for example, the stress field which is not subject to a Knudsen 
layer correction. 

Beyond these flaws, the major shortcoming of the work presented by Gibelli is that this 
author has avoided to compare his results with the outcomes obtained in some recent 
experimental studies. Indeed, the values of the slip coefficients determined in two inde- 
pendent series of experiments carried out for different gases |9(], [lO| show a very good 
agreement with those obtained in j3|, while there is a remarkable disagreement with the 



results presented in 



R 



concerning the second-order slip coefficient. 
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